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Here we present a problem related to the local Hamiltonian problem (identifying whether the 
ground state energy falls within one of two ranges) which is restricted to being translationally 
invariant. We prove that for problems with a fixed local dimension and 0(log(iV))-body local 
terms, or local dimension A^ and 2-body terms, there are instances of the problem which are QMA- 
complete. We discuss the implications for the computational complexity of finding ground states of 
these systems, and hence for any classical approximation techniques that one could apply including 
DMRG, Matrix Product States and MERA. One important example is a ID lattice of bosons with 
nearest-neighbor hopping at constant filling fraction i.e. a generalization of the Bose-Hubbard model. 



Introduction: The complexity class QMA is the quan- 
tum analogue of NP i.e. it is the class of problems for 
which there exists a polynomially sized quantum com- 
putation V that can verify the solutions to the problem, 
outputting 'yes' if the input is indeed a solution, and 'no' 
otherwise. Recent interest has focused around establish- 
ing complete problems for QMA, as these encapsulate the 
character of the class. This started with Kitaev's proof 
that 0(logiV)-local and 5-local Hamiltonian problems 
are QMA-complete where a Hamiltonian problem is 
defined to be establishing whether a particular Hamilto- 
nian has its smallest eigenvalue less than one threshold, 
a, or greater than another, 6, where b > a. Since this 
original proof, the structure of Hamiltonians which still 
solve QMA-complete problems has been reduced to 2- 
local Hamiltonians Q,and their spatial organisation in 
2D Q and ID lattices [J] has also been demonstrated. In 
the one-dimensional case, this was achieved by increasing 
the local dimension of spins from 2 to 12. Reductions 
of other, seemingly unrelated, problems have also been 
shown, hence proving their QMA-completeness, such as 
the A^-representability problem [j| . 

In this paper, we impose the translational invariance 
over a one-dimensional (periodic) chain on the Hamil- 
tonian, which can elucidate the trade-offs in resources 
required for quantum computation. One extreme that 
we already know is given by the aforementioned local 
Hamiltonian constructions [6|], where we have no time 
resolution in our control of the system, only spatial reso- 
lution. Another extreme is where we exchange the spatial 
control for temporal control 0], which is the culmina- 
tion of the study of global control schemes (see Q and 
references therein). The problem that we develop here, 
attempting to find a translationally invariant Hamilto- 
nian for computation, has neither spatial nor temporal 
control. Consequently, there is a cost elsewhere in the 
scheme - either in the local dimension of the spins used 
(~ poly(iV) where there are N qubits), or the range of 
the Hamiltonian terms (~ log(Af)). 

The ultimate aim of this study, however, is to exam- 
ine the computational complexity of finding the ground 
state energy of a translationally invariant local Hamil- 



tonian. It is typically argued that since there are ex- 
ponentially many parameters required to define a state, 
determining the ground state of an arbitrary Hamilto- 
nian is a hard problem. However, by restricting the sym- 
metry of a system, the effective number of parameters 
is reduced and thus one might hope to efficiently ap- 
proximate the ground state of, say, translationally invari- 
ant local Hamiltonians. This is a hugely important task 
in condensed matter physics since these are the Hamil- 
tonians that typically occur in nature, and is precisely 
the problem that is tackled by a range of approximation 
techniques such as DMRG/Matrix Product StatesjMPS) 
0, E3, Till- Quantum Monte Carlo simulations [l2j and 
MERA [13[ . There are good reasons to expect that effi- 
cient representations of these states exist [3, [l5| , and so 
it becomes a question of the computational difficulty of 
finding a good representation. One strikingly open ques- 
tion is "how accurately can the ground state of a one- 
dimensional translationally invariant Hamiltonian and 
the corresponding energy be approximated given poly- 
nomial time?", and it is hoped that one could discover 
techniques where it is possible to present certificates that 
guarantee a minimum level of accuracy. Some steps have 
been taken in understanding this problem, providing NP- 
hard instances of MPS [l6[ , and several classifications for 
the two-dimensional generalisation [Tt], EH • Developing 
QMA-complete instances of the translationally invariant 
local Hamiltonian problem will show that no approxima- 
tion strategy (that is applicable to these cases) can give a 

certificate without first resolving the question QMA=P. 
This is in contrast to the previous discussions [lg, E3, Ell 
which are tied to specific approximation strategies. 

Definitions: For completeness, we repeat the following 
definitions: 

Definition 1. (QMA) A promise problem L with 'yes' 
instances L yes and 'no ' instances L no is in QMA if there 
exists a quantum polynomial time verifier V and polyno- 
mial p such that for some e(M) — 2~ n ( M ) 

. Vx G L yes ,3 |0 Pr (V(\x) , |0) = 1) > 1 - e 

. VxGi no ,V|0 Pr(V(\x),\0) = l)<s 
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where \x\ = M, |£) e B® p{M) and Vn{V{\x) , |£» = 1) 
denotes the probability that V outputs 1 given \x) and 

Definition 2. Given a 2-local translationally invari- 
ant Hamiltonian on a one- dimensional chain of N D- 
dimensional systems, H = X)j=i with \\h\\ < 

poly(N), and two constants a < b, the promise problem 
that we examine has yes instances in which the smallest 
eigenvalue of H is at most a, and no instances in which 
it is larger than b, and we must decide which is the case. 
We call this problem TI2LH (Translationally Invariant 
2- Local Hamiltonian) . 

Computation from a Hamiltonian: The central part of 
previous proofs of QMA-completeness, as in this one, is 
the construction of Hamiltonians by showing how to im- 
plement an arbitrary quantum computation within the 
Hamiltonian (and hence can implement the verifier V). 
We shall assume that this computation acts on N qubits, 
of which the first M contain the string x to be veri- 
fied, and the other N — M are initialised in |0), and 
are used as ancillas in the computation. This compu- 
tation is typically achieved by introducing a 'clock' sys- 
tem, and a Hamiltonian that implements a particular 
gate as the clock increments. The Hamiltonian hops 
the system through the different clock states, in a way 
that maps directly to quantum state transfer along a spin 
chain [2(1 [2l[ . In the original proposal pj , the hop- 
ping terms between neighboring clock states all had the 
same strengths. Since a uniformly coupled chain does 
not achieve perfect transfer of a state through the sys- 
tem [Til [2l| (and hence perfect arrival of the computation 
in the final state is impossible), the computation can be 
extended by a number of identity operations, thus en- 
suring that the probability that the computation has got 
past the final step is high [2^|. It has since been ob- 
served @ by making an analogy with perfect state trans- 
fer schemes [2(1 Hl| that varying the coupling strengths 
can allow perfect arrival of the computation, and hence 
the addition of further steps is unnecessary. However, 
here we will continue to use the uniformly cou pled case 
due to the simplicity of the related eigenvectors [23[ , par- 
ticularly since for the proof of QMA-completeness, we do 
not actually perform the computation, and hence do not 
need perfect arrival of the result. 

Ultimately, to implement a computation composed of 
a discrete set of gates, we need to introduce a clock sys- 
tem. However, in a translationally invariant Hamilto- 
nian, we neither have time resolution nor spatial resolu- 
tion in which to encode this clock. Instead, we choose 
to hold it in the state of each spin. Each of these D- 
dimensional spins will have a number of states: on), 
and \x, y, z), where x G {0,1} is the qubit state which 
contains the computational component of the state, z G 
{1,2... N} denotes the position of the spin in the lat- 
tice and y G {0, 1 . . . R} is the time component, which we 
think of as 'cycle number'. The input state to our compu- 
tation will be a state of N+ 1 spins, \on) 0„ =1 \x n , n, 0) 



where the bits of the string x to be tested, and the 

\on) state indicates the initial position for a 'read/write- 
head' that our translationally invariant Hamiltonian will 
be able to scan backwards and forwards across the chain 
in discrete cycles (it will always scan across all N com- 
putational qubits). As it moves across, the cycle number 
is locally incremented by 1, and the required local uni- 
taries can be implemented due to the position label. The 
computing part of the Hamiltonian takes the form 

h% n » = -\<m) (o<®|0,l,0> <0,1,1L, 

-Un^m ® \m, n) (m - S rn , n\ i (g> |m, n') (m — 1 + 6 m , n'\. 

- U N. m ®\m,N- 1) (to, N - l\ { ® |to + 1, N) (m — 1,N 

~ u l% ® \ m + !> !> ( m - !) 1 U ® \m, 2} (m, 2\ d + h.c. 

where we are implicitly summing over cycle index to and 
position index n (n' = n + 1, 2S m = 1 + (—1)™). The 
unitaries U^ J m implement the step of the computation 
that is required to be acted at time to between qubits 
labelled n and n + 1, and only acts on the qubit space 
of sites i and j. The first term initialises the read/write 
head. The next term propagates the head to the right 
or left depending on the value of 8, and hence the par- 
ity of the cycle number. The final terms are responsible 
for reversing the direction of travel of the head when it 
reaches either extreme of the chain. For convenience, we 
fix the ground state energy of this Hamiltonian to zero 
by also adding the terms 

2V ® \m, n) (to, n\ t - V ® (|0, 1} (0, 1| + \R, 1) (R, 1|),. 

The final two terms mark a departure from the corre- 
spondence to transfer on a uniformly coupled chain, but 
makes the form of the lowest energy eigenvector of the 
Hamiltonian much simpler (see Eqn. ([T])). 

We now have a translationally invariant Hamiltonian 
Hcomp = where the terms \\h\\ < poly(N). 

However, it currently acts on a non-translationally in- 
variant input state. It is helpful to observe that the in- 
put states \rj n ), where \on) is located on qubit n can 
be combined to make the state translationally invariant, 
J2„ \ r l n )- After computation, the result would still be 
available because projective measurement would locate 
the position of \on). Alternatively, one could use the el- 
egant techniques described in [7[ to achieve computation 
on a translationally invariant input state. 

QMA-Completeness: Now that we have a translation- 
ally invariant Hamiltonian for performing computations, 
we turn our attention to proving the QMA-completeness 
of TI2LH. We do this by adding some extra terms to the 
Hamiltonian 

H = JlHi n p U t + J2{ctH f orm + H comp ) + R(N — 1)H output . 

Intuitively, the terms Hi npu t and Hf or m will verify that 
we have the desired input by adding an energy penalty 
for incorrect states. Hi npu t is particularly used to ensure 
the state of ancillas for the computation. We take these 
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to be positioned on qubits M + 1 to N, and must be 
in the |0) state at time 0. Thus, we require an energy 
penalty 

N 

h rn P ut = |l, ,n)(l,0,n| 4 . 

n=M+l 

The terms in Hf orrni while they could be applied just 
when the clock is at times 0, can equally well be ap- 
plied at all time because they act on parts of the state 
that don't change, such as the \on) and position states. 
As such, they are more conveniently considered in con- 
junction with H comp . The first term we apply is to en- 
sure that there is at least one |on) present in the system, 
— |on) (on\ i . However, we don't want more than one of 
these states, so we also need an energy cost for an |on) 
being on the right of a state in position n < N, that 
more than compensates for the reduced energy due to 
the | on) (on | term, 

>J 2 \x, m, n) (x, to, n| - <g |on) (on| 4+1 . 

n<N 

Similarly, we should ensure that neighboring position la- 
bels only increment by 1, 

22 \ x i m i n ) ( x i m ; n \i ® \ x > m,n + q) (x, m,n + q\ l+1 ■ 

Equally, H outpu t will test the final states, and add an 
energy penalty for a 'no' result, 

hr tput =11,^,1X1,^11,, 

where the computation outputs the accept/reject code 
(|0) / |1)) on the first qubit after R = poly(iV) cycles. 

We must now verify that the eigenstates of H satisfy 
the required properties of TI2LH. We choose to closely 
follow the proof presented in [2j , making use of their pro- 
jection lemma, 

Lemma 1. Let H = Hi + i? 2 be the sum of two Hamil- 
tonians operating on some Hilbert space H = S + S . 
The Hamiltonian H2 is such that S is a zero eigenspace 
and the eigenvectors in S have eigenvalue at least J > 
2||#i||. Then, 

m\s)-j^^<x(B)<xms). 

X(H) denotes the smallest eigenvalue of H . 

In particular, one can select, say J = 8\\Hi || 2 + 2\\Hi || 
in order to provide the bounds 

A(fli|s)-i<A(fO<A(fli| s ). 

Firstly, let's consider that there are 'yes' instances x 
with bits x n (and x n — if n > M). To test this solution, 
the input to the computation would have been 

N 

|?/o> = |on> \x n ,0,n) . 



We can consider the (N—1)R discrete steps of the compu- 
tation and write the state after each step as \r/ n ). Using 
these, we can write down a state 

and evaluate that 

(«| if input |n) = -1, (n| H comp |n) =0. 

Therefore, if the verifier accepts with probability > 1 — s, 
(n|ff|n)<-J 1+£ . 

Now we must prove that if there are no 'yes' instances, 
all eigenvalues are larger than some value b. One can 
precisely follow the arguments of 0. However, there 
is one small additional detail that we must first con- 
sider - the energy gap between the ground-state space 
of aHf orrn + H comp and the next excited state. Firstly 
we observe that due to Hf orm , the space splits into 
a series subspaces defined by the total number of |on) 
states present, and how badly ordered the position la- 
bels are. All of these subspaces have an offset of at least 
a from the original subspace of a single |on) state and 
correctly ordered position labels. Within this original 
subspace, we have the ground state subspace described 
by states |n) for all possible input values x, and cyclic 
permutations of the initial position of |on). Relative to 
the ground state space, all eigenvectors are separated 
by at least c/(R(N — l)) 2 for some positive constant 
c This is readily proved by performing a trans- 

formation to the state transfer model. Thus, by set- 
ting a > c/(R(N — l)) 2 , we can ensure that the mini- 
mum energy gap is at least c/(R(N — l)) 2 . We can now 
proceed with application of the projection lemma where 
H 2 = Ji{aHf OTm + H comp ) and i?i = H - H 2 . As in 0], 
we need to nest the application of this lemma, so that 
H = iii|scomp, H' 2 = J\Hi nput \s'"> m p and H[ = H — H' 2 
before finally arriving at the result that X(H) > 3/4 — e, 
thereby proving that solving the local Hamiltonian prob- 
lem implies the ability to solve all QMA problems. 

We can transform this Hamiltonian of N + 1 spins of 
dimension D = 2RN + 1 made up of two-body terms into 
a Hamiltonian acting on N log 2 (D—X)+l copies of 7-level 
spins and 2 log 2 (-D — l)-body terms. This transformation 
can be achieved by replacing each D-dimensional spin 
with a set of spins where one level is used as the |on) 
state, two levels encode a qubit state and the final three 
levels (3x2+1 = 7) denote whether that spin is being 
used to store a computational qubit, or a position or 
a cycle label, where an entire position and cycle label is 
achieved by combining log 2 (-D — 1) neighbors. We require 
N of these blocks, and a single extra spin to be initialised 
in the |on) state. The two-body terms from H comp need 
to act across two complete blocks of log 2 (D — 1) spins. 

Consequences and Conclusions: In summary, if there 
exist 'yes' instances to a QMA problem, the minimum 
eigenvalue of the corresponding H (implementing the ver- 
ifier V) is < — J 1 + e, whereas if there are no such in- 
stances, the minimum eigenvalue is at least —J\ + | — e. 
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Therefore, solving TI2LH allows you to solve any QMA 
problem. Note that our construction of the Hamiltonian 
is necessarily degenerate (24|, as seen by considering the 
case where there is a 'yes' solution, tested by \rj). There 
are N + 1 equivalent states corresponding to the different 
possible starting positions of \on), all of the same energy. 
A similar argument will hold for the case where there 
are no such solutions. Our construction also implicitly 
shows how one can take a system on ./V spins of dimen- 
sion d with /c-body local Hamiltonian terms and make it 
into a Hamiltonian acting on N spins of dimension Nd 
with /c-body local terms which is translationally invariant 
while preserving the structure of the eigenstates. 

This proves that any attempt to classically approxi- 
mate the ground state energy of a degenerate Hamilto- 
nian of N spins where either the Hamiltonian is composed 
of 0(log(A))-body (local) terms on spins of fixed dimen- 
sion, or 2-body terms acting on spins of poly(A) dimen- 
sion is QMA-hard (in N) and thus, assuming QMA^P, 
will, in the worst case, be intractable. Specifically, 
by polynomially increasing the computational resources 
available for determining the ground state energy, the ac- 
curacy, in the worst case, cannot be more than polynomi- 
ally increased. Falling within this class of Hamiltonians is 
the Bose-Hubbard model at constant filling fraction. At 
a filling fraction of 1, for example, if there are N lattice 
sites, then there are N bosons. We have to allow for the 
possibility that all N bosons could be found in a single 
lattice site, and thus each lattice site must be treated as 
an JV-level system. 

As we have proven, while there is a degenerate ground 
state, there are also excited states of these Hamiltoni- 
ans. Since there is an energy gap, properties such as 
exponentially decaying correlation functions [15j | hold, 
and techniques such as DMRG/MPS, and Vidal's re- 



cently proposed Multi-scale Entanglement Renormaliza- 
tion Ansatz (MERA) [l3[ provide an efficient representa- 
tion. We conclude that finding this representation must 
be a computationally hard task. This does not say, of 
course, that classical approximation strategies can't be 
useful in a range of physically important problems since 
the QMA-hardness of a problem only elucidates the worst 
case scenario. In contrast to previous attempts at similar 
problems [lj| [l?], Gil j the example Hamiltonian that we 
have produced is not tied to a specific algorithm that one 
might be trying to implement, but instead applies to the 
general question that we are trying to solve with such an 
algorithm, thereby eliminating the possibility of globally 
efficient algorithms that encapsulate the Hamiltonians in 
question. 

Evidently, an important question for the future is 
whether these results can be extended to the situation 
where both the local dimension and the Hamiltonian 
terms remain fixed as the system size grows, since this 
is the case which is most closely related to physical sys- 
tems. While it seems unlikely that our technique can be 
directly applied (as we need to use a clock), there are 
strong indications that this may well be the case. In the 
first case we analysed, where the local Hilbert space di- 
mension scales with poly (AT), but with two-body terms 
in the Hamiltonian, we know most of the structure of 
the ground state from Eqn. (fT]) if there are 'yes' solu- 
tions. By imposing this structure on the ansatz states of 
our approximation strategy, we are effectively reduced to 
minimising the energy over the local qubit states i.e. the 
QMA-hardness of this problem is effectively encoded in 
a nearest-neighbor Hamiltonian acting on qubits. 
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